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Abstract. We study the anomalous induced current of a vortex in a relativistic fluid via the 
chiral vortical effect, which is analogous to the anomalous current induced by a magnetic field 
via the chiral magnetic effect. We perform this analysis at weak and strong coupling. We 
discuss inequivalent implementations to the chemical potential for an anomalous symmetry. At 
strong coupling we use a holographic model with a pure gauge and mixed gauge-gravitational 
Chern-Simons term in the action. We discuss the holographic renormalization and show that 
the Chern-Simons terms do not induce new divergences. Strong and weak coupling results agree 
precisely. We also point out that the holographic calculation can be done without a singular 
gauge field configuration on the horizon of the black hole. 



1. Introduction 

Anomalies are responsible for the breakdown of a classical symmetry due to quantum effects. 
In vacuum the anomaly appears as the non-conservation of a classically conserved current in a 
triangle diagram with two additional currents. In four dimension two types of anomalies can be 
distinguished according to whether only spin one currents appear in the triangle [HE] or if also 
the energy-momentum tensor participates pJJ Hj . We will call the first type of anomalies simply 
chiral anomalies and the second type gravitational anomalies. In four dimension we should 
actually talk of mixed gauge-gravitational anomalies since triangle diagrams with only energy- 
momentum insertions are perfectly conserved (see e.g. [5]). In a basis of only left-handed fermions 
transforming under a symmetry generated by Ta the presence of chiral anomalies is detected by 
the non- vanishing of cIabc = ^ t (Ta{Tb,Tc}) whereas the presence of a gravitational anomaly 
is detected by the non- vanishing of bA = Ti(Ta)- 

Some studies showed that at finite temperature and density, anomalies give rise to new 
non-dissipative transport phenomena in the hydrodynamics of charged relativistic fluids. 
In particular magnetic fields in the fluid induce currents via the so-called chiral magnetic 
effect [SHE]- Later studies showed that a vortex in a fluid induces also a current parallel 
to the axial vorticity vector [91 [10]. On the basis of linear response theory, hydrodynamic 



transport coefficients can be extracted from the long-wavelength and low-frequency limits of 
some retarded Green functions. This leads to the so called Kubo formulas. For the chiral 
magnetic effect the Kubo formula has been derived in 112] . In [13] it was shown that the 
chiral vortical conductivity for charge and energy transport can be obtained respectively from 
the retarded Green functions 

° V = fc l™ o ^^ ea6c (J a T ob )U =0 , a^= kij^I^U, (1) 

c c a,fo c c a,fo 

where J % is the (anomalous) current and T %:> is the energy-momentum tensor (see also [141 [15] 
for details). 

In this manuscript we try to understand the effects anomalies have on the transport properties 
of relativistic fluids, both in the weak and strong coupling regimes, with special emphasis on the 
gravitational anomaly. Anomalies are very robust features of quantum field theories and do not 
depend on the details of the interactions. Therefore a non- interacting theory at weak coupling 
is sufficient for our purpose even without specifying to which gauge theory it corresponds to. 
By the same way a rather general model that implements the correct anomaly structure in the 
gauge-gravity setup is enough. Our approach for the latest case will therefore be a "bottom up" 
approach in which we simply add appropriate Chern-Simons terms that reproduce the relevant 
anomalies to the Einstein-Maxwell theory in five dimensions with negative cosmological constant. 



2. Weak Coupling 

In this section we compute the anomalous transport coefficients in a theory of free right-handed 
fermions tyf transforming under a global symmetry group G generated by matrices (T^y g . The 
chemical potential for the fermion tyf is given by \J = J^aQaI^a, where we write the Cartan 
generator Ha = Qa^ a m terms of its eigenvalues, the charges q A . We define the chemical 
potential through boundary conditions on the fermion fields around the thermal circle [16] 

*/( T ) = - e P^^f{j - p) } (2) 

with (3 = 1/T. The currents can be expressed in terms of Dirac fermions as 

N ■ N 

4= E^/^TV^, T 0i = ^^ / ( 7 °a i + 7 i 5V + ^, (3) 

/,<?=! f=l 

where we used the chiral projector V± = ^(1 db 75). The fermion propagator is 

= ^E A ^^W. A t (m;/,g)= 1 , (4) 

with i&f = iCj n + (1*, eft = (l,tq), q = E q = \q\ and uj n = 7rT(2n + 1) are the fermionic 
Matsubara frequencies. We can easily include left-handed fermions as well. 



2.1. Vortical conductivity 

The vortical conductivity is defined from the retarded correlation function of the current J\(x) 
and the energy momentum tensor or energy current T°i(x'), cf. Eq. ([3]), i.e. 

G V A (x - x') = \e ijn z8(t - f ) ([J A (x),TV(x')]) . (5) 
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Figure 1. 1 loop diagram contributing to the vortical conductivity Eq. ©. 



Going to Fourier space one gets the 1-loop contribution shown in Fig. [TJ The result for the zero 
frequency, zero momentum, vortical conductivity writes [17] 



The term involving the chemical potentials is induced by the chiral anomaly. More interesting 
is the term ~ T 2 with a coefficient that coincides with the gravitational anomaly coefficient [5]. 
This means that a non-zero value of this term have to be attributed to the presence of a 
gravitational anomaly. Left handed fermions contribute in the same way but with a relative 
minus sign. 

2.2. Magnetic conductivity 

The magnetic conductivity in the case of a vector and an axial U{\) symmetry was computed 
at weak coupling in [TT]. The corresponding Kubo formula involves the two point function of 
the current. Following the same method, we get the result for a general symmetry group 



No contribution proportional to the gravitational anomaly coefficient is found in this case. 
3. Chemical potentials for anomalous symmetries 

Before we go to how these results can be obtained at strong coupling from a gauge-gravity 
duality we stop for a moment and reflect on the formalism necessary to introduce the chemical 
potential. Note that we have been quite specific in how we introduced the chemical potentials 
in the weak coupling calculation in ([2]). Furthermore this is not the way chemical potentials are 
commonly discussed in textbooks. Rather than demanding twisted boundary conditions on the 
thermal circle it is far more common to consider a deformation of the Hamiltonian 



where Q is the charge in question. We can think of this as arising from the coupling of a (fiducial) 
gauge field to the current j M of the form / d^xA^j 11 and giving a vacuum expectation value to 
Aq = \i. With the fiducial gauge field we have gauge invariance now and we can remove of course 
the lj.Q coupling in the Hamiltonian by the gauge transformation A$ — > Aq + <9oX with \ = —Lit. 
Along the imaginary time direction t = —it this introduces of course just the twist in the 
boundary conditions on the fields in ([2]). As long as we have honest non-anomalous symmetries 
under consideration we have therefore two (gauge)-equivalent formalisms of how to introduce 
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Table 1. Two formalisms to chemical potential. 



Formalism 


Hamiltonian 


Boundary conditions 
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*( T ) = -*(r-/3) 
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the chemical potential summarized in table [T] [18J. One convenient point of view on formalism 
(B) is the following. In a real time Keldysh-Schwinger setup we demand some initial conditions 
at initial (real) time t = t^. These initial conditions are given by the boundary conditions in (B). 
From then on we do the (real) time development with the microscopic Hamiltonian H. This 
seems an especially suited approach to situations where the charge in question is not conserved 
by the real time dynamics. In the case of an anomalous symmetry we can start at t = t j with a 
state of certain charge but this charge does indeed decay over (real) time due to non-perturbative 
processes (instantons) or at finite temperture due to thermal sphaleron processes |19j . These 
processes are however supressed at large N and so can not be seen easily in the gauge-gravity 
correspondence. Taking this as excuse we simply ignore them, but keep them in the back of our 
head as motivation for favoring formalism (B) in the case of an anomalous symmetry, to which 
we come right now. 

Let us assume now that Q is an anomalous charge and start with our favoured formalism (B). 
We ask what happens if we do now the gauge transformation that would bring us to formalism 
(A). Since the symmetry is anomalous this means that the action transforms as 

S[A + d X ] = S[A] + J d 4 x X e^ pX (c^FpX + C 2 R a ^ U R? apX ) , (9) 

with the anomaly coefficients C\ and C2 depending on the chiral fermion content. It follows 
that formalisms (A) and (B) are physically inequivalent now, because of the anomaly. However, 
we would like to still come as close as possible to the formalism of (A) but in a form that is 
physically equivalent to the formalism (B). To achieve this we proceed by introducing a non- 
dynamical axion field 9(x) and the vertex 

Se [A, 9] = J d 4 x 9e^ A (ptF^F^ + C 2 R a ^ apX ) . (10) 

If we demand now that the "axion" transforms as 9 — > 9 — x under gauge transformations we 
see that the action 

S to t[A, 9] = S[A] + S e [A, 9] (11) 

is gauge invariant. Note that this does not mean that the theory is not anomalous now. We 
introduce it solely for the purpose to make clear how the action has to be modified such that two 
field configurations related by a gauge transformation are physically equivalent. In other words 
9 is a coupling and not a field. The gauge field configuration that corresponds to formalism (B) 
is simply Aq = 0. A gauge transformation with x = fJ>t on the gauge invariant action Stot makes 
clear that a physically equivalent theory is obtained by chosing the field configuration Aq = \x 
and the coupling 9 = —fit. If we define the current through the variation of the action with 
respect to the gauge field we get an additional contribution from Sq, 



j£=4C 1 e^ x d 1/ OF pX , 



(12) 



and evaluating this for = 



— fit we get the spatial current 



fg = 4C lf iB m , (13) 

(note that e 0l ^ k = — in Minkowski space). This is not the chiral magnetic effect! This is 
only the contribution to the current that comes from the new coupling that we are forced to 
introduce by going to formalism (A) from (B) in a (gauge)-equivalent way. The chiral magnetic 
and vortical effect are on the contrary non-trivial results of dynamical one-loop calculations as 
shown in the previous section. What is the Hamiltonian now based on the modified formalism 
(A)? We have to take of course the new coupling generated by the non-zero O. The Hamiltonian 
now is therefore^ 

H - fi (q + A J d?x (Cie^AidjAk + C 2 K°)j , (14) 

where K° is the zero component of the graviational Chern-Simons current = 
e^ x T% (d p T P aX + |r^r^ A ), fulfilling d^K* = \e^ x R a ^R? apX . 
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Figure 2. The three different contributions to the current current correlator relevant for the 
calculation of the chrial magnetic effect, (a) is the (UV-finite) finite temperature part, (b) is the 
part stemming from the (non-dynamical) axion vertex and (c) stems from the vacuum (T=0) 
triangle graph coupled to the background gauge field. 



Notice that in Formalism (A) we really have three contributions now to the current. One 
is formally a tree level contribution, namely the current (|13p . but there are also two different 
one-loop contributions. The first one comes from the UV-finite part of a finite temperature 
two point function. This is the part that is typically considered in the weak coupling Kubo 
formula calculations in I17j. But in formalism (A) we also have a (formally UV-divergent) 
contribution from the vacuum triangle diagram that couples to the gauge-field background! This 
has been first pointed out and emphasized in [21] . Graphically in formalism (A) we have to sum 
the graphs (a), (b) and (c) of figure [2j With = — fit and Aq = fi the contributions (b) and (c) 
cancel each other. 

We could also consider the contribution of the axion vertex to the energy-momentum tensor 
in the case when a gravitational anomaly is present. Since the Riemann tensor is however of 
second order in derivatives it is clear that the correspondent contribution to the energy current, 

This form of the Hamiltonian (without the gravitational part) is the one advocated in [20] , 



i.e. the T° l components of the energy-momentum tensor, will be of third order in derivatives^. 
In contrast the finite temperature one- loop graphs in Figs. [T] and [2] are first order in derivatives 
and contribute therefore to ordinary first order hydrodynamics! 

Our point of view is that the true anomalous transport effects are captured by the finite 
temperature graphs (a). These are also the only ones that contributes in formalism (B). We will 
apply these considerations now to the holographic calculation. 



4. Strong Coupling 

We present in this section the computation of the anomalous transport coefficients at strong 
coupling within a holographic model in five dimensions including terms which conveniently 
mimics the chiral and gauge- gravitational anomalies. 



4-1- Holographic Model 

Given an outward pointing normal vector n A oc 

an asymptotically AdS space with unit norm riAn A = 1, the induced metric takes the 
form fiAB = 9AB — n A n B- The action of our model is defined by 
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where Sgh is the usual Gibbons-Hawking boundary term and Da = /i^Vg is the covariant 
derivative on the four dimensional boundary. The second boundary term Scsk is needed 
if we want the model to reproduce the gravitational anomaly at general hypersurface. The 
most important fact about this action is the presence of two Chern-Simons terms. The one 
proportional to the parameter k is a pure gauge field CS term and the one proportional to A 
a mixed gauge-gravitational CS term. Both of them are diffeomorphism invariant, and they do 
depend however explicitly on the gauge connection Am- The bulk equations of motion are 

Gmn — AgMN = „FmlFn L 
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(19) 
(20) 



and they are gauge and diffeomorphism covariant. In Appendix A we discuss the holographic 
renormalization of the model within the Hamiltonian approach. This leads to the following 
counterterm of the action [TJ] 
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As a remarkable fact there is no contribution in 
the counterterm coming from the gauge-gravitational Chern-Simons term. This means that this 
term does not induce new divergences, and so the renormalization is not modified by it. 



2 This has recently been made explicit in [22] . 



4-2. Evaluation of Transport Coefficients 

The AdS/CFT dictionary tells us how to compute the retarded propagators [23} 124] . Since we 
are interested in the linear response limit, we split the metric and gauge field into a background 
part and a linear perturbation, 

9mn = 9mn + 6 h MN , A M = A$ +ea M ■ (22) 

Inserting these fluctuations-background fields in the action and expanding up to second order 
in e one can read the second order action which is needed to get the desired propagators [25] , 
If we construct a vector with the components of a,M an d Hmn an d Fourier transforming it, 
then it is possible to write the complete second order action on-shell as a boundary term 

5S ?l = I S^{&kAij& k J + &- k ®iM} ■ (23) 

J {2ir) a r-^oo 

From (j23|) it is possible to compute the holographic response functions by applying the 
prescription of [231 [2U [251 ES]. For a coupled system the holographic computation of the 
correlators consists in finding a maximal set of linearly independent solutions that satisfy infalling 
boundary conditions on the horizon and that source a single operator at the AdS boundary. To 
do so we construct a matrix of solutions F j(k, r) such that each of its columns corresponds to 
one of the independent solutions and normalize it to the unit matrix at the boundary. Finally 
using this decomposition one obtains the matrix of retarded Green functions 

Gu(k) = -2 lim (Aim(F M j(k, r))' + B u ) . (24) 

The system of equations f)19[> -(|20 |) admit the following exact background AdS Reissner- 
Nordstrom black-brane solution 

r 2 / \ T 2 MT 2 O 2 T 2 

d* 2 = ^(-/(^ + dx 2 ) +F A_d, 2 , f(r) = l-^ + ^, (25) 

A® = </>(f)dt = (a- ^ J dt, (26) 



where the horizon of the black hole is located at f = fn- We would like to draw attention to the 
boundary value of the gauge field a. We will make two choices, the choice (A) where a = \i and 
the choice (B) where a = 0. Note that with choice (A) the gauge field vanished at the horizon! 
At this point the reader will not be too surprised by our claim that the choice (A) corresponds 
to the formalism (A) of the previous section and choice (B) to the formalism (B). We know now 
however, that if we want (A) to be physically equivalent to (B) we also need to take into account 
the pure boundary action SeL4, 0]. In the holographic setup this is a pure boundary term, 
does not extend into the bulk and is therefore not dynamical. If we do not include this coupling 
it is known that the choices (A) and (B) give different results. In fact it has caused quite some 
confusion that the choice (A) without taking into account the ©-coupling gives a vanishing result 
for the chiral magnetic effect as shown in [8]. The alternative approach (B) causes however some 
uneasyness because now the gauge field necessarily is singular at the horizon (though not the field 
strength) , but it gives the correct result for the chiral magnetic effect with no need of including 
the pure boundary action SeLA, 0], as it was done in [21] . The nice feature of the modified 
approach (A) in holography is that we can now happily make the gauge field vanish at the 
horizon and still obtain the correct result for the anomalous conductivities on the boundary © 
The parameters M, Q and Hawking temperature of the RN black hole write 

M _4,Q 2 n _»rl T _ fl (2^M-3Q 2 ) 

M ~T 2 + 7^ ' Q -7T' T "4^ /(rH) - 2^ • (27) 



We just consider momentum fluctuations of the fields transverse to the perturbations of the 
momentum, i.e. we focus on the shear sector. Then one arrives at a system of four second order 
differential equations. The relevant physical boundary conditions on fields are: hf(0) = H a , 
B a (0) = B a ; where the 'tilde' parameters are the sources of the boundary operators. The second 
condition compatible with the ingoing one at the horizon is regularity for the gauge field and 
vanishing for the metric fluctuation [13] . 

After solving the system of differential equations perturbatively (see [H] for details), one 
can go back to the formula (|24p and compute the corresponding holographic Green functions. 
Finally, by using the Kubo formulas ([T]) one recovers the conductivities 

= -2^f = 4^' (28) 
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where we have included for completeness the chiral magnetic conductivity a 13 and the magnetic 
conductivity for energy current a 6 ' 13 . 

In a hydrodynamic framework we can summarize our findings in the constitutive relations 

T pu = (elP)uV + Prj pu + Q fl u u + Q*V, (31) 
J p = nu^ + N^, (32) 

with the first order in derivative terms 

N n = a ®B p + a v W, (33) 
Qi* = o^B* + <7 e ' v JF . (34) 



For simplicity of the expressions we have dropped here the usual dissipative terms related to shear 
and bulk viscosities or electric conductivity. The equilibrium quantities e,P,n are energy density, 
pressure and charge density and we defined the (covariant) mag netic field £> M = e^ px u v d p A x 
and vorticity vector ^ = e pupx u v d p u\. 

According to the discussion above, these values for the coefficients can be obtained, either 
by using the action of the holographic model and setting the deformation parameter a to zero 
(choice B), or considering the total action Stat[A, 0], Eq. (fTTjl . and setting a = \i and = — \it 
(choice A). In the latest case the axion term ()10p induces a contribution in matrices Au and Bu, 
see Eq. (|23p . The expression for a s is in perfect agreement with the literature and the one for a v 
shows the extra T 2 term already predicted in [17J and shown in Section [2] In fact the numerical 
coefficients coincide precisely with the ones obtained in weak coupling when specifying Eqs. ([6]) 
and ([7]) to the U(1)r group. This we take as a strong hint that the anomalous conductivities 
are indeed completely determined by the anomalies and are not renormalized beyond one loop. 



5. Discussion and conclusion 

In the presence of external sources for the energy momentum tensor and the currents, the 
anomaly is responsible for a non conservation law of the latter. This is conveniently expressed 
through [5] 



The axial and mixed gauge-gravitational anomaly coefficients are denned respectively by 



d A BC = ~[tT(T A {T B ,T c })R-tT(T A {T B ,T c }) L ], (36) 
b A = \x(T A ) R -tr(T A ) L , (37) 

where the subscripts R, L stand for the contributions of right-handed and left-handed fermions. 

We have computed the magnetic and vortical conductivity of a relativistic fluid at weak 
coupling and we find contributions that are proportional to the anomaly coefficients (j36[) and 
(j37|) . Non-zero values of these coefficients are a necessary and sufficient condition for the presence 
of anomalies [5] . Therefore the non- vanishing values of the transport coefficients ([6]) and ([7J) have 
to be attributed to the presence of chiral and gravitational anomalies. 

In order to perform the analysis at strong coupling via AdS/CFT methods, we have defined a 
holographic bottom up model that implements both anomalies in four dimensional field theory 
via gauge and mixed gauge-gravitational Chern-Simons terms. We have computed the anomalous 
magnetic and vortical conductivities from a charged black hole background and have found a 
non- vanishing vortical conductivity proportional to ~ T 2 . These terms are characteristic for the 
contribution of the gravitational anomaly and they even appear in an uncharged fluid. The T 2 
behavior had appeared already previously in neutrino physics [271 EE]- in |29| similar terms in the 
vortical conductivities have been argued for without any relation to the gravitational anomaly. 
However so far the effects of gravitational anomalies have not been taken into account in the 
purely hydrodynamic treatments, and therefore the T 2 terms appear simply as undetermined 
integration constants. The numerical values of the anomalous conductivities computed at strong 
coupling are in perfect agreement with weak coupling calculations, and this suggests the existence 
of a non-renormalization theorem including the contributions from the gravitational anomaly. 

So far we have computed the transport coefficients, and in particular their gravitational 
anomaly contributions, via Kubo formulas. It would be interesting to calculate directly 
the constitutive relations of the hydrodynamics of anomalous currents via the fluid/gravity 
correspondence within the model used in the present paper [91 [3Ql [3Tj . This study is currently 
in progress [32]. 
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Appendix A. Holographic Renormalization 

In order to have a deeper understanding of the holographic model, we will compute the 
renormalized action within the Hamiltonian approach, see e.g. |33|. 134] . Without loss of generality 
we choose a gauge in which A r = and the bulk metric writes 



ds z 



dr 2 + r jijdx l dx^ 



(A.l) 



Latin letters denote four dimensional (boundary) indices. The non vanishing Christoffel symbols 



are H, 



-K, 



Kj, where dot denotes differentiation respect r. Then one 
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can compute the off shell action in terms of transverse components of tensors. It is useful to 
divide the action up in three terms: 
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(A.4) 



The first one is the usual gravitational bulk and gauge terms with the usual Gibbons-Hawking 
term. Of particular concern is the last term in which contains explicitly the normal 

derivative of the extrinsic curvature Kij. Then the field equations will be generically of third 
order in r-derivatives. Having applications to holography in mind we can however impose the 
boundary condition that the metric has an asymptotically AdS expansion of the form 
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Using this expansion one can see that the last term in the action does not contribute in the 
limit r — > oo. Therefore the boundary action depends only on the boundary metric 7« but not 
on the derivative jij. 

The renormalization procedure follows from an expansion of the four dimensional quantities 
in eigenfunctions of the dilatation operator 
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For the extrinsic curvature and gauge fields, this expansion reads respectively 
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where the subindexes denote the corresponding eigenvalue, except for the cases 5dK(4) } = 
— 4-ff(4) *■ — 2K( 4 ) and <5d^4( 2 ) i = ~~ 2-A( 2 ) i — 2^( 2 ) j. Given the above expansion of the fields one 
has to solve the equations of motion in its Codazzi form, order by order in a recursive way. The 
derivative on r can be computed by using 



km 



2 1 d A xK l mllk 



km 



(A.9) 



By inserting in this equation the expansion of ifj given by Eq. ()A.7|) . one gets d/dr ~ 5d at the 
lowest order. Taking into account this, the computation of Kr \ *• is trivial, i.e. 
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By following the procedure and using the Codazzi form of equations of motion up to fourth 
order [13], one gets 
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In order to compute the counterterm for the on-shell action, besides the equations of motion 
an additional equation is needed. Following Ref. [34J, one can introduce a covariant variable 6 
and write the on-shell action as 
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The variable 9 admits also an expansion in eigenfunctions of <5d of the form 6 = 0r O ) + 9m\ + 
#(4) + #(4) loge _2r + • • •. Using the corresponding Codazzi equation for 9 [2], and the same 
procedure as above, one gets the result 
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'(2) 



(2 - d) ■ 

1 - 

4 



pipj _ P 2 _ U . .^ o) + 1 A ^ p _ ^p^ 



(A.14) 
(A.15) 



Finally the counterterm of the action can be read out from Eq. ()A.13j) by using K and 9 computed 
up to fourth order. The result is 



•>ct 



(d-1) 



1 + 



1 



(d-2)" 



1 



p»pj 
4(d-l) V J 4 



4^(0) ii^(o) y Jloge 



-2r 



(A.16) 



We have explicitly checked that the A dependence starts contributing at sixth order. So 
the gauge-gravitational Chern-Simons term does not induce new divergences, and there is no 
contribution in the counterterm coming from it. 
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